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Associated production of H ± and W T in high-energy e + e collisions 
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We study the associated production of the charged Higgs boson and W ± gauge boson in high 
energy e + e~ collisions in the Minimal Supersymmetric Standard Model (MSSM). This associated 
production, which first arises at the one loop level, offers the possibility of producing the charged 
Higgs boson at the e + e~ collider with mass more than half the center-of-mass energy, when the 
charged Higgs pair production is kinematically forbidden. We present analytic and numerical results 
for the cross section for e + e~ — > W + H~ in the full MSSM, taking into account the previously 
uncalculated contributions from supersymmetric (SUSY) particles. We find that the contributions 
of the SUSY particles enhance the cross section over most of SUSY parameter space, especially when 
the SUSY particles are light, ~ 200 GeV. With favorable SUSY parameters, at small tan /3, this 
process can yield more than ten W ± H^ 1 events for m H ± <J 350 GeV in 500 fb _1 at a 500 GeV e + e~ 
collider, or m H ± <, 600 GeV in 1000 fb" 1 at a 1000 GeV collider. 80% left-handed polarization of 
the e~ beam improves these reaches to m H ± <> 375 GeV and m H ± <J 670 GeV, respectively. 



PACS numbers: 12.60.Jv, 12.60.Fr, 14.80.Cp, 14.80.Ly 
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I. INTRODUCTION 



Discovery of the heavy Higgs bosons H°, A and -ff ± of 
the Minimal Supersymmetric Standard Model (MSSM) 1 
poses a special challenge at future colliders. Run II of 
the Fermilab Tevatron, now in progress, has a limited 
reach for the neutral heavy MSSM Higgs bosons. It can 
detect the CP-even and CP-odd neutral Higgs bosons 
H° and A if their masses are below about 150 GeV 
and tan /3 (the ratio of the two Higgs vacuum expecta- 
tion values) is large, so that the couplings of H° and 
A to b quarks are enhanced ^, ^, ||, || . For the charged 
Higgs boson H , no sensitivity is expected via direct pro- 
duction at the Tevatron unless QCD and SUSY effects 
conspire to enhance the cross section 0; however, i? ± 
can be discovered in top quark decays for m#± < m t 
and large tan/3 §, §. The CERN Large Hadron Col- 
lider (LHC) has a much greater reach for heavy MSSM 
Higgs boson discovery at moderate to large values of 
tan j3. H° and A can be discovered with decays to 
r pairs for tan/3 > 10 for a CP-odd Higgs mass of 
m A o = 250 GeV (tan/3 > 17 for m A o = 500 GeV) |, 0. 
The charged MSSM Higgs boson can be discovered 
in gb — > tH + with H + — > tv for virtually the same 
tan /3 values jl0|, [ll[] . The absence of a Higgs boson 
discovery at the CERN LEP-2 experiments implies that 
0.5 < tan/3 < 2.4 and m A o < 91.9 GeV are excluded 
at 95% confidence level |l2j . This leaves a wedge-shaped 
region of parameter space at moderate tan /3 in which the 
heavy MSSM Higgs bosons will not be discovered at the 



LHC. At a future high energy e + e~ linear collider (LC), 
the heavy Higgs bosons will be produced in pairs, if it is 
kinematically allowed. The dominant production modes 
are e + e~ — > H°A° and e + e~ — ► H + H~; for experimen- 
tal studies see Refs. jy| Q. These production modes 
are kinematically allowed only if m H o + m A o < ^J~s and 
2rriH± < y/s, respectively, where y/s is the center-of-mass 
energy of the LC. 2 At large m A o, m^o ~ m#o ~ mgi 
up to mass splittings of order m^/m^o, so that the 
pair-production modes are kinematically allowed only if 
771^40 < 0.5-y/s. In particular, the pair-production modes 
are limited to m A o < 250 GeV (to^o < 500 GeV) at a 
LC with Vs = 500 GeV (y/s = 1000 GeV). 

By contrast, the discovery of the light MSSM Higgs 
boson h , which typically has properties similar to those 
of the SM Higgs boson, is much more certain. At Run II 
of the Tevatron, discovery of h° at the 5a level is pos- 
sible for to/,0 < 120 GeV if 15 fb _1 of integrated lu- 
minosity can be collected and the detectors perform as 
expected [§. This 

covers a large fraction of the MSSM 
parameter space left unexcluded after LEP-2 At 
the LHC, discovery of h° is virtually guaranteed over all 
of the MSSM parameter space [^). Enough h° events 
are expected in a number of different production and de- 
cay modes to allow the measurement of various combina- 
tions of h° partial widths with precisions on the order of 
15% 0. At the LC, h° will be copiously produced. The 
promise of the LC for making precision measurements of 
the couplings of hP at the few-percent level has been well 
documented 0, ||, ||]. If the MSSM Higgs sector is 
not too far into the decoupling limit Jl9|, in which the 
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2 The pair production of scalars is p-wave suppressed near thresh- 
old, so in practice the Higgs mass reach in these modes is likely 
to be somewhat below the kinematic limit. 
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heavy Higgs bosons become increasingly heavy and the 
couplings of h° to SM particles approach their SM values 
(so that h° behaves like the SM Higgs boson), the LC 
measurements of hP couplings can be used to distinguish 
h° from the SM Higgs boson and to extract MSSM pa- 
rameters ^o). One finds > 2a deviations from the SM 
typically for m A o < 600 GeV |o|. 3 

While in the decoupling limit the couplings of h° 
become increasingly insensitive to the MSSM parame- 
ters, the couplings of the heavy Higgs bosons exhibit 
no such decoupling; rather, they are always sensitive to 
the MSSM parameters. Further, the measurements of 
the masses and couplings of the heavy Higgs bosons in 
addition to those of h° allow one to place orthogonal 
constraints on the parameters of the MSSM Higgs sec- 
tor. Thus, measurements of the properties of the heavy 
MSSM Higgs bosons are very valuable, especially in the 
decoupling limit. 

In this paper we consider the production of one of the 
heavy Higgs bosons alone or in association with lighter 
SM particles at the LC. While the cross sections for such 
production modes are typically very small, they offer the 
possibility of extending the reach of the LC to higher 
values of m A a H a H ± > 0.5^. Single heavy Higgs boson 
production has been studied in the context of the MSSM 
or a general two Higgs doublet model (2HDM) in a num- 
ber of processes. The following final states have been 
considered in e + e _ collisions: ZH° and h°A° 
and e+e-H° (see^e.^Ref. p||); bbH°, bbA\r 



and tbH+ (23 
and W+H- 




9lH; 



24J2 

HT Wa 153]. If the e + e LC is converted 
into a photon collider through Compton backscattering 
of intense laser beams, single heavy Higgs bosons can be 
produced via 77 collisions. In particular, H° and A can 
be produced in the s-channel pj, |35| , |36| ; the final states 
t~vH + and tbH + |57j and W^H^ |3§[1 are also accessi- 
ble. Finally, in e~7 collisions one can produce vH~ |39[ . 

In this paper we compute the cross section for e + e~ — > 
W + H~ in the full MSSM. We present analytic formu- 
lae for the matrix elements and numerical results for the 
cross section as a function of the MSSM parameters. We 
confirm the calculation of e + e~ — > W + H~ in the 11011- 
supersymmetric 2HDM performed in Refs. |31 ' ! ) ! !l 1 



We find that the contributions of the SUSY particles en- 
hance the cross section over most of SUSY parameter 
space, especially when the SUSY particles are light, with 



3 A subset of the h° couplings can also be measured with similar 
precision at a LC operating as a 77 collider, leading again to 
> 2a deviations from the SM typically for m A o ^ 600 GeV [ p]J . 

4 A significant contribution to e + e — — > W + H~ comes from the 
loop-induced W + H~ Z and W + H~^ vertices. Ihe contributions 
to these vertices from top/bottom quark loops pn] and gauge and 
Higgs boson loops |4l[ have been computed in the 2HDM. The 
full MSSM contributions to these vertices have been computed 
in Ref. Jf^ |; the contributions from top/bottom squark loops are 
also given in Ref. pol. The loop-induced W+H~Z vertex has 
also been computed numerically in the full MSSM in Ref. Ejl . 



masses of order 200 GeV. The largest contributions to 
the cross section from the SUSY sector come from di- 
agrams involving charginos/neutralinos and top/bottom 
squarks. We also find that left-handed polarization of 
the e~ beams leads to an order 50% enhancement of the 
cross section. Throughout this paper we present cross 
sections for the single process e + e~ — > W + H~; these 
numbers should be doubled to find the combined cross 
sections for this process plus its charge conjugate. 

This paper is organized as follows. In Sec. [Il] we in- 
troduce our form alism for the matrix elements and cross 
section. In Sec. Ill we display the relevant SUSY di- 
agrams and review the renormalization procedure. In 



Sec. IV we present our numerical results. In Sec. [v| we 
review the literature on other single heavy Higgs boson 
production processes and compare their reach to that of 
e + e~ — > W + H~ . Sec. VI is reserved for our conclusions. 



The matrix elements for the 2HDM and SUSY diagrams 
are collected in the appendix. 



II. FORMALISM 

Following the notation of Ref. [ 3l] , the matrix element 
M. for e + e~ — > W + H~ can be decomposed into six in- 
dependent matrix elements Ai and their corresponding 
coefficients M.^. 



i=l 

The six matrix elements are defined as, 

A! 



(1) 



A 2 
A 3 
A, 
A 5 
A 



V(j>2) fa 1 "t, 75 u (Pl)(Pl ' e *( fc l)) 



v{pi) fa 
v(Pz) fa 
v(pz) fa 



2 

1-75 
2 

1+75 

2 

1-75 



u(pi)(pi ■ e*Oi)) 
u(Pi)(P2 ■ e*Oi)) 
u(Pi)(P2 ■ e*(&i)), 



(2) 



where e* is the polarization vector of the W + boson, pi 
and P2 are the incoming momenta of the initial e~ and 
e + , respectively, and fci is the outgoing momentum of the 
W+. 

For convenience we define two additional matrix ele- 
ments, A7 and As- 



A 7 
A 8 



is 



a ^Mp2h»^^u( Pl y a k 10 k 2l , (3) 



where k 2 — P\ + p 2 — k\ is the outgoing momentum of 
the H" and £ q ^t is the totally antisymmetric tensor, 
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with e 0123 = 1. At and A s can be expressed in terms of 
Ai~Aq as follows: 



At = Ai(t - u)/2 + A 3 - A 5 
A s = A 2 (u-t)/2 - Ai + A 6 , 



(4) 



where t = (pi — ki) 2 and u = (pi — k 2 ) 2 . 

From the form of these matrix elements, one immedi- 



J 



ately sees that the cross sections for like polarizations of 
e + and e~ are zero: 



d<J ( e R e R -> W+H-) _ da(e\el -> W+H~ 



e?cos 6 



dcosO 



= 0. (5) 



For unlike polarizations of e + and e , the cross sections 
are given in terms of the six basic matrix elements by: 



da(e+el -> W+H~) 



H 



2s\M: 



,2 W± ffl ff 



tu) 



dcos9 327rs 



Am? w 



[(m 2 w - t) 2 \M A \ 2 + (m 2 w - u) 2 \M e \ 2 



A(m 2 w - t)Re[M 2 M* 4 ] + 4(m 2 w - u)Re[M 2 M* 6 ] + 2(tu - m 2 w s - m 2 w m 2 H± )Re[M 4 M* 6 ] + 4|M 2 | 2 }] , (6) 



and 



da(ete R ^W+H-) 



2 S \M 1 \ 2 - { ^ ±m2w 



±m 2 w 



dcosO 32irs 
k(m 2 w - t)Re[MiM* 3 ] + A(m 2 w - u)Re[MiM* 5 ] 



{(m 2 w - t) 2 \M 3 \ 2 + (m 2 w - U ) 2 |.M 5 | 2 



+ 2(tu - m 2 w s - m 2 v m 2 H± )Re[M 3 M* 5 ] + 4|A^i| 2 | 
I 



(7) 



Here s = (pi + p 2 ) 2 and k = 2\ki\/^/s is given by: 



K 2 = 



1 



(m H i 



m w ) 



(m H ± - m w f 



(8) 



The unpolarized cross section is obtained by averaging 
over the four possible initial combinations of e + e~ polar- 
izations. 



III. RENORMALIZATION 



The diagrams contributing to e + e — » W + H 



the 

3 1H. We show the addi- 
Thcre are also W~ H~ 



2HDM are shown in Fig. 
tional SUSY diagrams in Fig. 
and G~ H~ mixing self-energy diagrams that involve SM 
fermions, gauge and Higgs bosons (Fig. ||), and SUSY 
particles (Fig. ||). These mixing self-energies, together 
with the counterterms, contribute to e + e~ — > W + H~ 
via the diagrams in Fig. |^. 

We neglect all diagrams that are proportional to the 
electron Yukawa coupling. We also neglect the diagrams 
shown in Fig. ^: 

(a) These two diagrams are proportional to m e either 
through the electron Yukawa coupling or through 
the factor of m e obtained via the equation of mo- 
tion of the incoming electron. In particular, the 
vector boson-Higgs mixing is proportional to = 
Pip, + P 2 fj,, which gives m e when acting on the 
V^ej^e vertex. 

(b) This diagram is zero because the W~H~ mixing 
self-energy is proportional to fci„, and ki • e* = 0. 









(12) 



FIG. 1: Feynman diagrams for the 2HDM contributions to 
e + e~ — > W + H~. Here S, Sij,k denote Higgs and Goldstone 
bosons, V = 7, Z, and V' = Z,W ± . 



(c) We will set the renormalized tadpoles to zero below, 
so that this diagram does not contribute. (Note 
that the A and G° tadpoles are zero automatically 
due to CP conservation.) 

(d) This diagram is purely real and is canceled by the 
G~ H~ mixing countcrterm, as discussed below. 

We evaluate the one-loop process e + e~ — » W + H~ in 
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e+ H- e+ H- e~ Xi H~ 




e + (10) c- (n) W* 



FIG. 2: Feynman diagrams for the SUSY contributions to 
e + e~ — * W + H~ . Here Xi,j,fc denote charginos and neutrali- 
nos, Si,j,k denote squarks and sleptons, and V = 7, Z. 



b Si Si 




(1) (2) (3) 



b S, S 




(4) (5) (6) 



FIG. 3: Feynman diagrams for the 2HDM contributions to 
W~ H~ and G~ H~ mixing self-energies. Here Si = h°,H° 
and Sj = H ± ,G ± . 



Xi s t 




FIG. 4: Feynman diagrams for the SUSY contributions to 
W~ H~ and G~ H~ mixing self-energies. Here xt,j denote 
charginos and neutralinos and Si,j denote squarks and slep- 
tons. 




FIG. 5: Feynman diagrams for the contributions to e + e~ — > 
W + H~ through the counterterms and the W~ H~ and 
G~H~ mixing self-energies. The "X" in diagram 1 denotes 
the counterterm for the 'yW + H~ and ZW + H~ vertex, and 
the blob in diagrams 2-4 denotes the renormalized W~H~ 
and G~ H~ mixing. 



H- e+ H- 




II- 



FIG. 6: Contributions to e e — > W H that we neglect. 

the 't Hooft-Feynman gauge using dimensional regular- 
ization. Using dimensional reduction yields the same 
result for this process. Diagrams that contain closed 
loops of chiral fermions may give rise to inconsistencies 
in the treatment of 7 5 in naive dimensional regulariza- 
tion ^5) . We have checked that this does not happen 
in the one-loop e + e~ — > W + H~ diagrams by comput- 
ing the terms involving 7 s in the diagrams with closed 
fermion loops using the 't Hooft-Veltman rules ff"6|| . 
It was shown in Ref. [Q that this prescription yields a 
consistent formulation of dimensional regularization even 
if 7 5 couplings are present. (For a pedagogical discussion 
see Ref. [^J.) We find that in this process the 't Hooft- 
Vcltman rules yield the same result as naive dimensional 
regularization, so that no inconsistencies arise. 

We follow the on-shcll renormalization formalism de- 
veloped by Dabelstein Q for the MSSM gauge and 
Higgs sectors. The one-particle-irreducible (1PI) tadpole 
diagrams for h° and H° are given by — iTh and —iTn, 
respectively. Including the tadpole counterterms, the 
renormalized tadpoles are given by: 

f h =T h +5t h , f H = T H + 5t H . (9) 

The 1PI two-point function for W + H + mixing is de- 
fined as —ik^Yi W + H + (k 2 ), where k is the incoming mo- 
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mentum of the and H + is outgoing. 5 The conjugate 
two-point function, with W~ incoming and H~ outgoing, 
is given by +ik^Ti W + H + (k 2 ), where again k is the incom- 
ing momentum of the W~ . The renormalized two-point 
function for W + H + mixing is obtained by adding the 
counterterm: 

^w+H+{k 2 ) = T, w+H +(k 2 ) 

( 5v 5v \ 

—mw sin (3 cos (3 I 5Zr — 5Zh, 1 ) , (10) 

V vi v 2 J 

where 5Zu t and 5Zh 2 renormalize the two Higgs doublet 
fields H\ and H 2 and Svi and 5v 2 renormalize the two 
Higgs vacuum expectation values Q : 

Hi -> z]lfHi\ Vi -> Z]lf{vi - 5vi), (11) 

and Zhi = 1 + 5Z^ i . 

The 1PI two-point function for G + H + mixing is de- 
fined as +iT, G + H + (fc 2 ), where k is the incoming mo- 
mentum of the G + , and H + is outgoing. The conju- 
gate two-point function, with G~ incoming and H ~ out- 
going, is the same. The renormalized two-point func- 
tion for G + H + mixing, T,q+ h + (fc 2 ), is fixed in terms of 
Yi W+H+ (fc 2 ) due to the Slavnov- Taylor identity (see, e.g., 
Refs. @ |0) for details): 

k 2 t w+H + (fc 2 ) - mw^ G +H+ (fc 2 ) = 0. (12) 

Finally, there are two vertex counterterms for 
^pW^H^ and Z^W^H" (all particles outgoing). 
These are given by — iemw sin /3 cos (3 Scg^ and 
+iesw'mz sin (3 cos (3 5c g^ u , respectively. Here 5c = 
(5Zn 1 — 5Zh 2 — 5vi/vi + 5V2/V2) and sw denotes the 
sine of the weak mixing angle. Note that 5c is fixed in 
terms of the W + H + mixing counterterm, Eq. |Io| . 

We need impose only the following two renormalization 
conditions. 6 First, the renormalized tadpoles in Eq. || 
are set to zero. Second, the real part of the renormalized 
W + H + mixing is set to zero when H + is on mass shell: 

Re± w+H+ {m 2 H± ) = 0. (13) 

This fixes the following combination of counterterms: 

(5v 5v \ 
5Zh 1 — 5Zh 2 1 I 
vi v 2 J 

= ReZ w+H+ (m 2 H± ), (14) 

which appears in the jW + H~ and ZW + H~ vertex 
counterterms. Applying the Slavnov- Taylor identity 



5 We use the convention X> M = c* M + ieA^. Rcf. |31| uses the 
opposite convention, so our matrix elements shouladiffcr from 
theirs by a sign. 

6 Because the process e + e — — > W + H~ is zero at tree level, the 
renormalization procedure is greatly simplified and many renor- 
malization conditions, such as that for tan/3, need not be im- 
posed. 



(Eq. [12]) at fc 2 = rn 2 H± , this condition also fixes 
Re (m 2 H± ) = 0. In addition, since the countert- 

erms are purely real, we have from Eq. p"2| : 

TYl 2 

Im S g +h+ {m 2 H± ) = — ^Im T, w+H + (m 2 H ±), (15) 
mw 

so that diagrams 4-6 of Fig. ||, diagrams 3 and 4 of Fig. |] 
and diagram (d) of Fig. |^ need not be calculated. As a 
check of our calculation we have verified Eq. |l5| explicitly, 
including the full MSSM contributions. 

Explicit results for H w +h+(j> 2 ) are given in the ap- 
pendix. The real part of Y, w + H + (m 2 ^± ) fixes the 
-fW + H~ and ZW + H~ counterterms. The imaginary 
part of Yjw+h+ ( m #-± ) remains in diagrams 1 and 2 of 
Fig. ||, since it is not canceled by the W + H + mixing 
counterterm. This imaginary part also determines dia- 
gram 3 of Fig. H via Eq. [l^. 

We have checked explicitly that all the divergences in 
the 1PI diagrams that contribute to e + e~ — > W + H~ are 
canceled by the counterterms. 

IV. NUMERICAL RESULTS 

In this section we examine the behavior of the cross 
section for e + e~ — * W + H~ for various choices of SUSY 
parameters and evaluate the regions of parameter space 
in which the cross section is large enough to be observed 
for m H ± > y/s/2. We assume data samples of 500 fb _1 
at y/s = 500 GeV and 1000 fb" 1 at y/s = 1000 GeV. We 
choose an optimistic standard of detectability to be ten 
H ± production events in the LC data sample. Adding 
together the cross sections for W + H~ and W~H + pro- 
duction, our standard corresponds to a cross section for 
W+H- production of 0.01 fb at -y/s = 500 GeV and 
0.005 fb at V« = 1000 GeV. We assume that the e+ 
beams are unpolarized. We compare the cross sections 
obtained with unpolarized e~ beams with those obtained 
with 80% left- or right-polarized e~ beams. We do not 
make any attempt to consider backgrounds or apply cuts; 
this is beyond the scope of our present analysis. 

Unlike the case of the non-supersymmetric 2HDM, in 
which the top/bottom quark loops give by far the largest 
contribution to the cross section, in the full MSSM the 
fermionic loops involving charginos/neutralinos and the 
bosonic loops involving stops/sbottoms also give contri- 
butions of similar size. Although the stop/sbottom loops 
are enhanced by the large H~t^b* L coupling (which is 
proportional to the top quark Yukawa coupling), these 
diagrams are suppressed by higher powers of the super- 
particle masses than the fermionic loops. 

Diagram 1 in Fig. |^ and diagrams 1 and 3 in Fig. ^ de- 
couple in the limit of heavy gauginos and Higgsinos, while 
diagrams 3-10 in Fig. || decouple in the limit of either 
heavy sleptons or of heavy gauginos/Higgsinos. Squarks 
and sleptons contribute to diagram 2 in Fig. I and 2 and 
4 in Fig. |], where stops/sbottoms give the largest contri- 
bution because of the large top quark Yukawa coupling. 
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In our numerical analysis, we use a common sfermion 
mass scale MgusY = 200 GeV for all the squarks and 
sleptons. Msusy is the soft SUSY-breaking mass param- 
eter that enters the diagonal elements of the squark and 
slepton mass matrices. We also include the additional D- 
term contributions to squark and slepton masses, which is 
crucial for the Slavnov- Taylor identity (Eq. [l5]) to hold. 
We study two different choices for the trilinear A cou- 
plings: (I) A t = A b = and (II) A t = A b = 200 GeV. 
Along with the fj, parameter and tan/3, A t determines the 
left-right mixing in the stop sector, which plays an im- 
portant role for relatively light Msusy- In addition, the 
H parameter determines the Higgsino masses. In most of 
our analysis we choose \i = 500 GeV; we also consider 
fx = 100 GeV when examining the tan/3 dependence of 
the cross section. We fix the U(l) and SU(2) gaugino 
mass parameters to be 2 Mi = M2 = 200 GeV. Working 
consistently at the one-loop level, we use the tree-level re- 
lations for the Higgs masses and mixing angles in terms of 
m,A° and tan/3. We have verified numerically that using 
the radiatively corrected values for the CP-even Higgs 
masses and mixing angle does not change our numerical 
results in any significant way. 

In what follows, we have taken all the SUSY break- 
ing masses for the squarks and sleptons to be the same 
for simplicity. With Msusy as low as 200 GeV, the ra- 
diatively corrected mass of the lightest CP-even MSSM 
Higgs boson h° lies below the current experimental Higgs 
search bound. However, it is of course possible to choose 
Msusy for the stop/sbottom sector to be large and/or 
to impose large left-right mixing in the stop sector so 
that the mass of h° is increased above the experimental 
bound, while still keeping the sleptons and first two gen- 
erations of squarks relatively light so that they give large 
contributions to the e + e~ — > M /± _ff =F cross section. In 
this case, the numerical results presented below will of 
course change slightly, but the general conclusions from 
our analysis will remain true. A more detailed analy- 
sis of the constraints due to the h° mass bound will be 
presented elsewhere pl[ . 

Figures [j] and |^ show the dependence of the e + e~ — > 
W + H~ cross section on the charged Higgs mass for 
tan/3 = 2.5 and 7, respectively. In both cases we plot 
cross sections for yfs = 500 GeV (left) and -s/s = 1000 
GeV (right). Solid lines are the contributions from the 
non-SUSY 2HDM with the Higgs sector constrained by 
the MSSM mass and coupling relations. 7 The dotted 
(dashed) lines show the cross sections in the full MSSM, 



7 The pure 2HDM contributions to e+e — — > W + H~ have been 
studied in Refs. ]3l| Our numerical results for the 2HDM 

are in close agreement with those of Refs. fell B3] , after taking 
into account a factor of 1/4 from the average over the initial e + e — 
polarizations that was omitted in Ref. [pl|, and noting that the 
results in Ref. Jj^ | are the average over the spin states ej~e~j 
and e^e - ^; for unpolarized beams the cross sections in Ref. |32| 
should be divided by two. 



including the contributions from all the superparticles, 
for A t = A b = (A t = A b = 200 GeV). We compare 
the cross sections with 80% left-handed e~ polarization, 
no polarization, and 80% right-handed e~ polarization, 
which are denoted in each plot by the same type of lines, 
from top to bottom. Left-handed e~ polarization always 
gives a larger cross section. The additional SUSY con- 
tributions generally enhance the cross section. In certain 
cases, for example -^=1000 GeV and tan/3 = 7 (Fig. ||), 
the SUSY contributions increase the cross section by al- 
most an order of magnitude. The cross sections decline 
as rn.jj± increases; however, reasonable cross sections can 
be obtained for m H ± > \fs/2, especially for small tan/3. 

The ten-event discovery reach in m#± is shown in Ta- 
ble U for a 500 GeV LC (requiring a cross section of 0.01 
fb) and a 1000 GeV LC (requiring a cross section of 0.005 
fb). Using an 80% left-polarized electron beam generally 
increases the reach by at least 20 GeV compared to the 
unpolarized case. The increase in the reach in the MSSM 
compared to the 2HDM is larger at larger tan j3 or higher 
LC center-of-mass energy. In the following discussion, we 
will consider the results for an 80% left-polarized electron 
beam. For tan/3 = 2.5, a charged Higgs with mass up to 
375 GeV (667 GeV) could be detected at a 500 GeV (1000 
GeV) LC. For favorable SUSY parameters, the SUSY 
contributions can increase the reach by about 20 GeV at 
v/s = 500 GeV and by more than 100 GeV at ^ = 1000 
GeV, compared to the non-SUSY 2HDM. At tan/3 = 7, 
the cross section is generally too small to be observed at 
the ten-event level; however, a reach in m H ± up to about 
580 GeV is still possible at a 1000 GeV machine. The de- 
pendence of the cross section on A t for low Msusy = 200 
GeV is due to the effects of left-right mixing in the stop 
sector. For larger values of Msusy, the squark contribu- 
tions become less important and the dependence on A t 
becomes much less significant. The cross section is also 
sensitive to the value of the /i parameter. For example, 
lowering fi to 100 GeV increases the ten-event discovery 
reach in m#± to 390 GeV at tan (3 — 2.5 or 300 GeV for 
tan/3 = 7 at a 500 GeV LC, almost independent of A t 
(see Fig. |). 

Figure ^ shows the tan/3 dependence of the e + e~ — > 
W + H~ cross section for 80% left-polarized electrons. We 
plot two different values for /i, 100 and 500 GeV, to show 
the [i dependence as well as the dependence on A t for 
small Msusy = 200 GeV. For n = 500 GeV, values of 
tan j3 above 20 lead to a negative mass-squared for the 
lighter sbottom and are forbidden. From Fig. || we can 
see that while the cross section in the 2HDM (solid line) 
falls rapidly with increasing tan /3, the cross section in the 
MSSM experiences a much milder drop as tan increases. 
Thus the reach in tan/3 in the MSSM is larger than in 
the 2HDM. At a 500 GeV machine, for m H ± = 350 GeV, 
the ten-event discovery reach is tan /3 < 5.8 in the MSSM 
with favorable parameters as shown in Fig. ^|, compared 
to tan/3 < 2.5 in the 2HDM. At a 1000 GeV machine, for 
ro#± = 600 GeV, the ten-event discovery reach is tan /3 <, 
5 in the MSSM, while in the 2HDM even tan/3 - 2 is 
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FIG. 7: The e + e" -> W + H~ cross section as a function of m H ± for tan/3 = 2.5, at s/s = 500 GeV (left) and 1000 GeV (right). 
The trilinear couplings are chosen as At — At — (dotted lines) and 200 GeV (dashed lines). The rest of the SUSY parameters 
are chosen to be Msusy = 200 GeV, 2Mi = M2 = 200 GeV, and /j, — 500 GeV. The solid lines show the cross section in the 
non-SUSY 2HDM (with MSSM relations for the Higgs sector). In each plot, the lines from top to bottom are the cross sections 
for 80% left-polarized, unpolarized, and 80% right-polarized electrons. 




TABLE I: Ten-event discovery reach in m H ± (in GeV) at a LC with an 80% left-polarized electron beam, for ^J~s = 500 GeV 
(C = 500 fb _1 ) and 1000 GeV (C = 1000 fb _1 ). The corresponding reach with an unpolarized electron beam is given in 
parentheses. The SUSY parameters are as specified in the caption of Fig. W. 





tan/3 = 2.5 


tan/3 = 7 




500 GeV 


1000 GeV 


500 GeV 


1000 GeV 


2HDM 

MSSM, At = At = 
MSSM, At = A b = 200 GeV 


352 (327) 
362 (344) 
375 (347) 


535 (<500) 
540 (512) 
667 (600) 


< 250 

< 250 

< 250 


< 500 

< 500 
581 (517) 
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FIG. 9: The e + e~ — > W + H~ cross section as a function of tan/3 for 80% left-polarized electrons, with y/s = 500 GeV, 
m H ± — 350 GeV (left) and y/s — 1000 GeV, m H ± = 600 GeV (right). Four different choices for \i and A t = ^6 are shown: 
fi = 500 GeV, A t = A b = (dotted line); ^ = 500 GeV, A t = A b = 200 GeV (dashed line); [i = 100 GeV, A t = At = 
(dash-dotted line); and /i = 100 GeV, At = A b = 200 GeV ("+" line). The remaining SUSY parameters are Msusy = 200 
GeV and 2Mi = M 2 = 200 GeV. The solid lines show the cross section in the non-SUSY 2HDM (with MSSM relations for the 
Higgs sector). 



not detectable. In Table || we show the maximum ten- 
event discovery reach in tan/3 when m H ± is just above 
y/s/2. For a 500 GeV machine, we find that a reach up 
to tan/3 < 9 is possible for m H ± ~ 250 GeV. For a 1000 
GeV machine, while in most cases we find a reach up to 
tan j3 < 5 for m#± ~ 500 GeV, in certain cases a wide 
range of tan/3 values can be explored. 

For the relatively low values of Msusy, Mi and M^ 
used in our analysis, the dependence on is complicated. 
We find an enhancement of the cross section over the 
2HDM for most values of fi. However, there are values 
of \i that lead to a large suppression in the cross section. 
For example, at yfs = 500 GeV with a left-polarized elec- 
tron beam, m H ± = 350 GeV, tan/3 =2.5, A t = A b = 200 
GeV, and fi around 800 GeV, there is a large cancel- 
lation between the 2HDM and SUSY matrix elements, 
leading to a suppression of the cross section. Further- 
more, changing the relative signs of Mi, Mi and \x can 
alter the SUSY contribution to the cross section. 



V. COMPARISON WITH OTHER CHANNELS 

Single heavy Higgs boson production has been stud- 
ied before via a number of different processes, which we 
summarize here. Because detailed experimental studies 
of almost all of these processes are unavailable, we again 
choose an optimistic standard of detectability to be 10 
heavy Higgs boson production events in the LC data sam- 
ple. We assume data samples of 500 fb -1 at yfs = 500 
GeV and 1000 ib" 1 at -y/s = 1000 GeV. For neutral Higgs 
boson production, this 10-event standard corresponds to 



a cross section of 0.02 fb at ^/i = 500 GeV (0.01 fb at 
yfs = 1000 GeV). For charged Higgs boson production, 
we add together the cross sections for H + and H~ pro- 
duction before applying the 10-event standard. In what 
follows we assume that the e + and e~ beams are unpo- 
larized, and adapt the cross sections presented in the lit- 
erature accordingly. We consider only tan /3 values above 
the LEP lower bound of 2.4 (I2) and heavy Higgs masses 
above y/s/2. 

A. e + e~ collisions 

The 10-event reach for various single heavy Higgs bo- 
son production modes in 500 GeV and 1000 GeV e + e~ 
collisions is shown in Figs. [l(] and [ll], respectively. 

The heavy CP-even neutral Higgs boson H° can be 
produced via Higgsstrahlung, e + e~ — » ZH°; H^-boson 
fusion, e + e~ — * vDH ; and Z-boson fusion, e + e~ — > 
e + e~H°. The cross sections for these processes are sup- 
pressed relative to the corresponding SM Higgs produc- 
tion cross sections by cos 2 (/3 — a), where a is the mixing 
angle that diagonalizes the CP-even neutral Higgs bo- 
son mass-squared matrix. In the decoupling limit |l9|| , 
cos 2 (P — a) oc m% / m\ , so these cross sections decrease 
rapidly as m^o increases. At a 500 GeV LC, ZH° pro- 
duction yields > 10 events well beyond m H o ~ 0.5y/s, out 
to m H o < 380 GeV and for low to moderate tan (3 < 25 
(see Fig. |l0]). The contour shown in Fig. [h] was produced 
using the program FeynHiggsXS p^] , which includes the 
full Feynman-diagrammatic corrections to this process 
at one-loop as well as two-loop corrections to the Higgs 
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TABLE II: Ten-event discovery reach in tan/3 for m H ± ~ \fs/2 at a LC with an 80% left-polarized electron beam, for y/s = 500 
GeV (£. = 500 fb _1 ) and 1000 GeV (€. = 1000 fb _1 ). The corresponding reach with an unpolarized electron beam is given in 
parentheses. The SUSY parameters are as specified in the caption of Fig. H. 
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8.7 (5.9) 
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3.8 (2.9) 
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FIG. 10: Ten-event contours for single heavy Higgs boson 
production in unpolarized 500 GeV e + e~ collisions, for 500 
fb _1 of integrated luminosity. Solid lines show H° produc- 
tion via ZH° (22) and bbH° g (the bbH° line is a 20-event 
contour: see text). Long-dashed lines show A production via 
h°A° g and ZA° E|. Short-dashed lines show H+ pro- 
duction via tbH + and tvH + also shown are our results 
for W + H~ production in the 2HDM (lower curve) and full 
MSSM with light superpartners (upper curve). On the i-axis 
we plot the mass of the relevant heavy Higgs boson. See text 
for details. 



masses and mixing. 8 At a 1000 GeV LC, the cross sec- 
tion for ZH° production is too small to be of use. The 
cross section for Higgs boson production via VF-boson fu- 
sion typically becomes more important as y/s increases 
(see, e.g., Ref. fllcf ). While the cross section for vvH° at 
y/s = 500 GeV is too small to be of use, at yfs = 1000 
GeV it yields > 10 events for rn#o < 550 GeV and 
tan/3 < 4.5 (see Fig. O). 9 Finally, the Z-boson fusion 



The FeynHiggsXS |22| input parameters were chosen as 
^SUSY = 1000 GeV for all three generations, M 2 = 200 GeV, 
M s = 500 GeV, fj, = -200 GeV, and A t = A b = 2M SUSY to 
yield the maximal mixing scenario. Mi was fixed by the GUT 
relation to M2. 

These results for e+e — — > uPH° are ba sed on the c* c ff approxima- 
tion from the program subhpoldm |52| , with the same MSSM in- 
put parameters as used in the FeynHiggsXS computation above. 
A full Feynman-diagrammatic calculation of the radiative cor- 
rections to this process is not yet available. 
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FIG. 11: As in Fig. |T| but for 1000 GeV e+e" collisions and 
1000 fb _1 of integrated luminosity. The solid line shows ui>H° 
production. The long-dashed line shows h°A° production p2| . 
Short-dashed lines show H + production via ibH + and ri>H + 
p5[ ; also shown are our results for W + H~ production in the 
full MSSM with light superpartners. 



cross section is always about a factor of ten smaller than 
the H^-boson fusion cross section, so we do not expect it 
to be useful here. 

In addition to the kincmatically limited production 
mode e + e~ — > H°A°, the heavy CP-odd neutral Higgs 
boson can be produced in association with the light h , 
via e + e~ — > h°A°. As above, the cross section for this 
production mode is suppressed by cos 2 (/3 — a). Using 
FeynHiggsXS p2| with input parameters as given above, 
we find that at a 500 GeV LC h°A° production yields 
> 10 events for m^o < 350 GeV and low to moderate 
tan 13 < 14 (see Fig. 0); at 1000 GeV the 10-event reach 
is m A o < 560 GeV and tan^ < 3.5 (see Fig. 0. 

The heavy MSSM Higgs bosons can also be produced 
in association with pairs of third-generation fermions: 
e+e- -> bbH°, bbA° , T-DH+, and ibH+ . The cross sec- 
tions for the first three of these processes are strongly 
enhanced at large tan /3, and the fourth is enhanced at 
both large and small tan f3. 10 At a 500 GeV LC, the pro- 



10 The first two processes are useful for measuring tan /3 at large 
tan /3 and low m A a . 
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cess e + e~ — > is especially promising: it yields > 20 

events for m H o < 360 GeV and large tan/3 > 10 [^3|, 
this 20 event contour is shown in Fig. [lOl 11 The process 
e + e~ — > 66A° yields less than 20 events for m^o > 250 
GeV [§| H]. At a 1000 GeV LC, neither bbH° nor bbA° 
production yields more than 10 events for m#o ^40 > 500 
GeV §|. 

At a 500 GeV LC, both T -pH+ and tbH+ production 
yield > 10 events at large tan (3 ~ 40 for m#± < 270 
GeV, with t~DH + production having a slightly larger 
cross_ section g, || |6) (see Fig. 0). At a 1000 GeV 
LC, tbH + production is more promising, due to the larger 
phase space available; it yields > 10 events at large 
tan/3 ~ 40 for m#± < 550 GeV, while t~vH + produc- 
tion gives a reach of only m#± < 520 GeV {23L |25|] (see 
Fig. @). 12 

Finally we consider single heavy MSSM Higgs produc- 
tion modes that are zero at the tree level but arise at one 
loop. The process e + e _ — > ZA° has been calculated in 
both the 2HDM @ and in the full MSSM ||. 13 In the 
2HDM the cross section is too small to be of use at ei- 
ther Vs = 500 or 1000 GeV |§7). In the full MSSM, light 
charginos and sleptons with masses of order 200 GeV 
can enhance the cross section by an order of magnitude 
or more; this leads to > 10 events at a 500 GeV LC for 
m A o < 270 GeV and tan/3 < 5 || (see Fig. [lC]). This 
enhancement disappears once the chargino and slepton 
masses exceed about 500 GeV. For y/s = 1000 GeV, even 
with light SUSY particles the cross section is too small to 



be of interest here. The similar process e 



7 A was 



considered in Refs. |£7J, [29, |30j; however, even including 
the enhancement of the cross section from the contribu- 
tions of light SUSY particles, this process yields too few 
events to be of use at either yfs = 500 or 1000 GeV. 

Comparing the process e + e~ — > W + H~ to the various 
processes described above, we see that W + H~ produc- 
tion is the only channel in e + e~ collisions analyzed to 
date that yields > 10 events containing charged Higgs 
bosons at low tan/3 values for > \/~s/2- In the 

2HDM with MSSM relations imposed upon the Higgs 
sector, W + H~ yields > 10 events for m H ± < 325 GeV 
and tan/3 < 4.5 at ^fs = 500 GeV (Fig. [h]), while for 
\fs = 1000 GeV the cross section in the 2HDM is too 
small to be of use. If the MSSM constraint on the Higgs 



11 Results were given in Refs. [^, |24j for tree-level e+e — — > bbH° , 
bbA° cross sections down to 0.1 fb, which yields 50 signal 
events, using mj = 4.25 GeV. We take into account the dom- 
inant QCD corrections jj^ by scaling these cross sections by 
{m b (m H+ )/ml° lc ) 2 ~ 0.4 for ~ 250 - 350 GeV, so that 
the cross section of 0.1 fb in Refs. fea, K4] corresponds to roughly 
20 signal events after QCD corrections are included. 

12 The QCD corrections to e+e — — > F6_ff+ were recently computed 
in Ref. [55) and were found to reduce the cross section at large 
tan /3; in particular, for t/s = 500 GeV and tan j3 = 40, the QCD 
corrections reduce the 10-event reach by about 8 GeV [B5| . 

13 A significant contribution to e+e - — > ZA° comes from the loop- 
induced A°ZZ vertex, computed in Ref. |5q|. 



sector is relaxed, then the e + e~ — » W + H~ cross sec- 
tion can be enhanced in some regions of 2HDM param- 
eter space due to large triple-Higgs couplings Un- 
der the requirement that the Higgs self-couplings remain 
perturbative, this enhancement can increase the reach at 
Vs = 500 GeV up to tan /3 < 6.5 for m H ± < 280 GeV 
Q |§. Including the full MSSM contributions, 14 the 
reach increases to m H ± < 345 GeV and tan j3 < 5 at 
= 500 GeV (Fig. 0), and m H ± £ 580 GeV and 
tan/3 < 3.5 at ^Js — 1000 GeV (the latter reach is com- 
parable to that of h°A° production; see Fig. pi] ). 



B. 77 and e 7 collisions 

If the e + e _ LC is converted into a photon collider 
through Compton backscattering of intense laser beams, 
the neutral heavy Higgs bosons H° and A can be singly 
produced in the s-channel through their loop-induced 
couplings to photon pairs. This process appears to be 
very promising for detecting H° and A° with masses 
above \/s/2 and moderate tan/3 values between 2.5 and 
10 [[34], |35|, p6[ . In particular j_a recent realistic simulation 
of signal and backgrounds |34| showed that a 630 GeV 
e + e~ LC running in 77 mode for three years would al- 
low H°, A detection over a large fraction of the LHC 
wedge region (in which the heavy MSSM Higgs bosons 
would not be discovered at the LHC) for m^o up to the 
photon-photon energy limit of ~ 500 GeV. At a 1000 
GeV LC, the mass reach is likely to be above 600 GeV 
II- 

The cross sections for production of r~ vH + and tbH + 
in 77 collisions |5?J are expected to be larger than the 
corresponding cross sections in e + e~ collisions at large 
tan/3. In particular, if 1000 fb -1 of e + e~ luminosity at 
y/s ee = 1000 GeV is devoted to 77 collisions, t~DH + 
production yields > 10 events for m H ± up to almost 700 
GeV at large tan /3 ~ 30 0. 

Production of W + H~ in 77 collisions also occurs at 
the one-loop level. The cross section for this process has 
been computed in the non-supersymmetric 2HDM and 
the full MSSM in Ref. @. If 1000 fb" 1 of e+e~ luminos- 
ity at 1000 GeV is devoted to 77 collisions, > 10 events 
would be produced in the 2HDM for m H ± < 570 GeV at 
tan/5 = 2; the cross section falls with increasing tan/3. 
In the full MSSM with light superparticles the cross sec- 
tion is enhanced, and this mode provides > 10 events for 



The MSSM parameters were chosen here as MgygY = 1000 GeV 
for the third generation squarks, MguSY = 200 GeV for the 
rest of the squarks and sleptons, M 2 = 2Mi = 200 GeV, A t = 
Af, = 2000 GeV (i.e., the maximal mixing scenario), and fi = 200 
GeV for = 500 GeV and fi = 500 GeV for y^s = 1000 GeV. 
The tree-level relations were used to determine the MSSM Higgs 
masses and couplings; we checked numerically that using the 
radiatively-corrected MSSM Higgs masses and couplings does not 
change our results significantly fell]. 
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m H ± above 600 GeV at tan/3 = 2 (m H ± < 520 GeV 
at tan/3 = 6) 38 . Assuming comparable luminosities, 
77 — > W + H~ is competitive with e + e~ — > W + H~ in 
the MSSM. This implies that spending some of the e + e~ 
luminosity on running in photon collider mode would not 
significantly impact the number of W + H~ signal events 
collected. 

Finally, if the LC is run in e~7 mode, the process 
e~7 — > vH~ is possible. The cross section for this pro- 
cess has been computed in the 2HDM |39| and found to 
be relatively independent of m#± all the way up to the 
kinematic threshold for parent e + e _ center-of-mass ener- 
gies of 500 or 1000 GeV. Unfortunately, with the typical 
expected e~7 luminosity of 100 fb _1 [ p7) , the cross sec- 
tion for this process in the 2HDM is too small to be of 
interest |]39| . This process could become promising in the 
MSSM if its cross section is enhanced by the contribu- 
tions of light superpartners |5jJ| , or if the e~7 luminosity 
is increased. 



VI. CONCLUSIONS 

We studied the single charged Higgs production pro- 
cess e+e- -> W ± HT at a 500 GeV and a 1000 GeV 
linear collider, including the complete MSSM contribu- 
tions at one loop. At small tan/3, the reach in to#± 
from this process extends above the kinematic threshold 
(m H ± ~ y/s/2) for e+e~ — > H + H~ pair production. We 
found that light SUSY particles with masses of order 200 
GeV enhance the cross section in most of the parameter 
space. At small tan /3, an increase in the reach in m#± by 
about 20 GeV in the full MSSM compared to that in the 
2HDM is possible at a 500 GeV machine, and by more 
than 100 GeV at a 1000 GeV machine. In particular, with 
favorable SUSY parameters and tan (3 = 2.5, this process 
yields more than 10 W^H^ events for m H ± < 350 GeV 
in 500 ft)" 1 at a 500 GeV collider, or m H ± < 600 GeV 
in 1000 fb _1 at a 1000 GeV collider, assuming unpolar- 
ized beams. At large tan/3, an order of magnitude en- 
hancement of the cross section compared to the 2HDM is 
possible, so that the cross section suffers only a moder- 
ate decrease with increasing tan/3. The reach in tan /3 is 
therefore enhanced; with favorable SUSY parameters and 
m H ± ~ y/s/2, this process yields more than 10 W ± H T 
events for tan/3 < 6 at a 500 GeV collider. We also 
found that an 80% left-polarized electron beam enhances 
the cross section by about 50% and improves the discov- 
ery reach by at least an additional 25 GeV; in particular, 
the reach at tan/3 = 2.5 improves to m H ± < 375 GeV at 
a 500 GeV collider, or m H ± < 670 GeV at a 1000 GeV 
collider, and the reach in tan/3 improves to ~ 8.5 at a 
500 GeV collider. 

The process e+e~ -> W + H~ in the MSSM with light 
superpartners is very promising compared to the other 
single charged Higgs boson production modes that have 



been considered to date. At the 



collider, W+H~ 



tbH + production at high tan/3. At an e~7 collider, the 
cross section for vH~ production in the 2HDM is too 
small to be seen with a typical integrated luminosity of 
100 fb _1 in the high-energy part of the 7 spectrum; how- 
ever, this process could become promising in the MSSM 
if its cross section is enhanced by the contributions of 
light superpartners, or if the e~j luminosity is increased. 
Finally, if 1000 fb" 1 of e+e" luminosity at 1000 GeV is 
devoted to 77 collisions, the process 77 — > W + H~ in the 
MSSM will be competitive to e+e~ -> W + H~. 

A few additional processes that have not yet been 
computed may be promising for single heavy Higgs bo- 
son production at an e + e~ collider. The behavior of 
the SM Higgs production cross sections leads us to ex- 
pect that weak boson fusion processes will naturally have 
larger cross sections than Higgsstrahlung-type processes 
at y/s ~ 1000 GeV. The process e + e _ — > De~H + may 
thus be promising J58|. The process e + e~ — > vvA° may 
also be promising; however, preliminary results for the 
top/bottom quark loop contributions in the 2HDM in- 
dicate that this process is too small to be observed for 
m A a > y/s/2 Q. This process could become observable 
if the addition of the full 2HDM and MSSM contributions 
enhances the cross section sufficiently p8l. 
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APPENDIX A: NOTATION AND CONVENTIONS 

For couplings and Feynman rules we follow the con- 
ventions of Ref . JEp . We follow the notation of Ref . |3l[] 
for the one-loop integrals. The one-point integral is: 



16tt 2 



A(m 2 



d D q 



1 



(2n) D (q 2 - m 2 ) 



(Al) 



where D is the number of dimensions. The two-point 
integrals are: 



167T' 



{B Q ,k^B x }{k 2 ,mimi) = 



(A2) 



production at low tan /3 is complementary to r vH + and 



d D q 



{W} 



{2n) D {q 2 -ml){{q + k) 2 -m\) 



■2\ ■ 
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The three-point integrals are: 



16?r 2 
d D q 



{C ,C^C^} = 



(A3) 



{l,q»,q»q»} 



(2ir) D (q 2 - m 2 )((q + ^i) 2 - m 2 )((<? + fe) 2 - to 2 ) ' 



where the tensor integrals are decomposed in terms of 
scalar components as 



= k^Ci + k%C 2 



+(k^ + k^)C 12 . 



(A4) 



The arguments of the scalar three-point integrals are 
(kf, (k 2 — fci) 2 , k\, to 2 ,, mf, to 2 ,). The four-point integrals 
are: 



167T* 



{D^D^D^} 



{l,q^q^} 



(2n)D ( q 2 - m 2 )((q + h) 2 - m 2 )((q + k 2 ) 2 - to 2 )((<z + fc 3 ) 



to 2 ) : 



r 



(A5) 



where the tensor integrals are decomposed in terms of 
scalar components as 



k%D 2 



KD* 



g^Doo + kfkZDn + k^k u 2 D 22 + k^k^D. 
+ (fcf k% + kgkfiDia + {k»k» + k%k»)D 13 



(A6) 



The arguments of the scalar four-point inte- 



grals are (k 2 , (k 2 - 

fc3) 2 ,TO§,TO?,m!,TO§). 



ki) 2 ,(k 3 - k 2 y,kiMX^ 



APPENDIX B: 2HDM CONTRIBUTIONS 

Explicit results for the matrix elements for e + e~ — > 
W+H- in the 2HDM were given in Ref. 0. We have 
calculated all the 2HDM diagrams independently and 
give the formulae here for completeness. After correct- 
ing a few typographical errors ]60| in the formulae of 
Ref. pH , we find agreement with all of their results ex- 
cept for Figs. ||.2 and |[3, which differ by an overall sign. 

The quark triangles that appear in Figs, [jj.l and [l].2 



were also computed in Ref. p0| and agree with our re- 
sults. (Note that Ref. jl0| uses the convention £0123 = 1, 
which is opposite to our convention.) 

For convenience, we list here our notation for the gauge 
and Yukawa coupling coefficients used in this section. 

The photon coupling coefficients to leptons/quarks are: 



4 L 



4 R 



(Bl) 



where the electric charges are e„ = 0, e e = — 1, e u = 2/3, 
and e e i = —1/3. 

The Z boson coupling coefficients to leptons/quarks 
are: 



9z 



' 1 - (~T 3 + efSw)/swcw, 9 f z R = (efs^/swcw, 



(B2) 

where T3 = 1/2 for v, u and T3 = —1/2 for e, d. 

For the W boson coupling to leptons/quarks we define: 



9w 



= -1/V2. 



sw- 



im 



The and G !± coupling coefficients to top/bottom 
quarks are (the couplings to leptons and first and second 
generation quarks are small and taken to be zero): 



Vntb 


Vntb 


Vatb 


Vatb 


rat cot (3/V2mwsw 


nib tan/3/ \[2mwsw 


—mb/V^mwsw 


mt / V^mw sw 



The coefficients for couplings of three Higgs bosons are: 



9h 


-AG+ 


—imw/2sw 


9h- 


-h°H+ 


—mwspa/sw — mz cos 2(3 sin(/3 + a)/2swcw 


9h- 


H ( >H+ 


—m\vCf3 a /sw + mz cos2/3cos(/3 + a)/2swCw 


9h- 


-h°G+ 


mwCf3a/2sw — mz sin 2/3 sin(/3 + a)/2swcw 


9h- 


H°G+ 


—mwS[3a/2sw + mz sin 2/3 cos(/3 + a) /2swcw 


9g- 


-h°G+ 


mz cos2/3sin(/3 + a) /2s\\rcw 


9g- 


H°G+ 


—mz cos 2(3 cos(/3 + a) /2swcw 
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where sp a = sin(/3 — a) and cp a — cos(/3 — a), and we have imposed the MSSM relations on the H + H~h°, H + H~H° 
couplings. 

We now list our results for the 2HDM diagrams. 
Fig. |l: 



M 



2a 2 N c g w 



{[-(mtg^ymb + ™ b (g v L - 9v R )ymb)Bo + 2g v L (m t y% tb + m b y% tb )C 00 



+ {~9v L {mty L H tb + mbVmbWw + ( m t9v L Vmb ~ ™b<?yVL) fc i • ht)C x 
+(9v L (m t ymb + mty^h ■ k 2 - (m t g v L y^ tb ~ m b g v R y§ tb )m 2 H± )C 2 

+ (-™hv L VHtb - m 2 t m b (g v L - g v R )y§ tb + m t m 2 g v R y L mb + m t gfy L mb k x ■ k 2 )C Q ][g e v R Ui + g e v L A 2 ] 
+ [-2gf(m t y L mb + m b y R tb )(C 12 + C 22 ) - gf{im t y L mb + m b y R tb )C 2 - {m t gfy L mb - m b g* R y R tb )G\ 
-m t gir L y^ tb Co}lgv R (A 3 + A 5 ) + g^ L (A A + As)] 

+ [mtgv L ym b Co + {m t 9v L ymb + ^bSy^iJ^i + 9v L ( m tymb + ™bV% tb )C 2 }(g e v R A 7 + g e v L A 8 )} , 

with the arguments for the integral functions as -B(s, m 2 , m b ), C(m^, s, mj I± , mf , m b , m b ). 
Fig. 0.2: Similar to Fig. | .1 with the exchange of 



(B4) 



mt*->m b , 9v L ^gv L , 



dR uR 

9v ^ 9v ! 



ymb VHtbl A*, A s terms flip sign. 



(B5) 



Fig. 0.3: 

M = 



a 



-gvwwgwwsgnws 



x { [-B Q - (s + m 2 - m 2 w )C - (a + m 2 w - m 2 H± )d + (s - m 2 w + m 2 H± )C 2 + C 00 ] + g e v L A 2 ) 

+(4Ci + C 2 - C 12 - C 22 ) [g c v R (A 3 + As) + g e v L (A 4 + A e )] } , (B6) 

with the arguments for the integral functions as B(s, m|p, m|^), C(m^, s,m%±,mg,m 2 v ,m 2 v ). The couplings are 
given in the following table: 



s 


gjww 


gzww 


9HWS 


gwws 


h° 


-l 


— Cw/ Sw 


cp a /2sw 


m W S{3a/sw 


H° 


-l 


—c w / sw 


— Sf3 a /2s W 


mwcpa/sw 



Fig. 0.4: 



M = 



2a 2 



TO' 



-9zSiSj9HwSigwwSj {Coo(gy A\ + g^A 2 ) 



(2C + 2d + 3C 2 + Ci 2 + C 22 ) [g c v R (A 3 + A 5 ) + g v L {A A + A)] } 



(B7) 



with the arguments for the integral functions as C{m 2 v , s, rn 2 H ± , i 
following table: 



5 . , m 2 s . ) . The couplings are given in the 



Si Sj 


gzSiSj 


9HWSi 


gwws. 


A h° 
A H° 


-icp a /2swcw 

iSf3a/2s W CW 


—i/2sw 
—i/2sw 


mwspa/sw 
m w cp a /s w 



Fig. 0.5: 



M 



-gvs z s 3 gHs i s k gws j s k {Coa(g v R Ai + gy L A 2 ) 
(C 2 + Cia + C 22 ) [g e v R (A 3 + A 5 ) + g e v L (A 4 + A e )] } , (B8) 

g fc ,m|.,TOg.). The couplings are given in the 



with the arguments for the integral functions as C(m 2 v ,s,mij±,rir- < .m 
following table: 



14 



s t 


Si 


s k 




gzSiSj 




9WSjS k 


A 


h° 


G ± 





-icp a /2swcw 


9h 


-AG+ 


— Sj3a/2sw 


A 


H° 


G ± 





isp a /2s w cw 


9h 


-AG+ 


— Cf3 a /2sw 


h° 


A a 


H± 





iCf3a/2s W CW 


9h- 


h°H+ 


i/2sw 


H° 


A {} 


H± 





-isp a /2s w cw 


9h- 


H°H+ 


i/2s w 


h° 


G° 


G ± 





iSj3a/2s W Cw 


9h- 


-h°G+ 


i/2sw 


H° 


G° 


G ± 





icp a /2s w c w 


9h- 


H°G+ 


i/2s\y 


G ± 


G ± 


h° 


1 


cos(29w) /2swCw 


9h- 


-h"G+ 


sp a /2sw 


G ± 


G ± 


H° 


1 


cos(29w) /2swcw 


9h- 


H°G+ 


Cf3 a /2SW 


H± H± 


h° 


1 


cos(29w) /2swcw 


9h- 


h°H+ 


Cf3 a /2sw 




H° 


1 


cos(29w) /2swcw 


9h- 


H°H+ 


— Sf3a/2sw 



Fig. 0.6: 

2 

M = T9vvs t gwv>s ] gHs t s ] Co{gv R Ai +gy L M), (B9) 

s — my 

with the arguments for the integral functions as C(m^, s,mj I± ,mg.,m 2 / , ,m|.). The couplings are given in the 
following table: 





Si 


V 


9jVSi 


gzv'Si 


gnSiSj 


gwv'Sj 


h° 


G ± 


z 





mzspa/swcw 


gH-h ( >G+ 


—mzsw 


H° 


G ± 


z 





m z cp a / s w cw 


9h—h°g+ 


-m z s w 


G ± 


h° 


w ± 


mw 


—mzsw 


9H-h°G+ 


mwspa/sw 


G ± 


H° 


w ± 


mw 


-mzsw 


9h-h°g+ 


m w cf} a /s w 



Fig. 0.7: 

2a 2 

M = 2'9vvs i gHV'S :j gws i s j {-Caoig^Ai + gv L A 2 ) 

s — m v 

+(-C 2 + C12 + C 22 ) [g v R (A 3 + A 5 ) + gf{M + A)] } , (BIO) 

with the arguments for the integral functions as C(m\ v , s,m 2 H ±,m 2 s .,m 2 Si ,rn v ,). The couplings are given in the 
following table: 



Si 


s s 


V 


9jvs z 


gzv'Si 


gnv'Sj 


gwsiSj 


h° 


H± 


Z 





m z sf3 a /s w cw 


cos(29w)/2swcw 


— Cf3 a /2sw 


H° 


H± 


z 





m z cp a /swcw 


cos(29w)/2swcw 


Sj3 a j2sw 


G ± 


h° 


w ± 


mw 


—mzsw 


C[3a/2SW 


S[3 a /2SW 


G± 


H° 


w ± 


mw 


-m z s w 


—sp a /2sw 


cp a /2sw 



Fig. 0.8: 

KA a " 

•A4 — n-gzzsgHwsgzww 

s — m z 

{ [B + (2m 2 H± + 3m 2 w - 2s)C - (s - m 2 H± - 3rr&)Ci - (s - m 2 w - 3m 2 H± )C 2 - C 0Q ] (g v R A 1 + g e v L A 2 ) 
+ (2Co-2C 1 +3C 2 + C 12 + C 22 )[g v R (A 3 +A 5 )+gy L (A4+A 6 )]}, (Bll) 

I 



with the arguments for the integral functions as 
B(s,m 2 y 1 m'g), C(m 2 v ,s,m 2 H± ,m 2 v ,m 2 z ,m 2 s ). The cou- 
plings are given in the following table: 



s 


gzzs 


gnws 


gzww 


h° 


mzspaj 's w cw 

m Z C{3a/ S W CW 


Cf3 a /2sw 
— Sf3 a /2sw 


—c-wl sw 
—cwl sw 
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Fig. 0.9: 

M = 



<x 9vwhs9w ws 

S — TO?/ 



with the arguments for the integral functions as 
B (m^ , m s , rriyy) . The couplings are given in the fol- 
lowing table: 



5 


g^Hws 


9ZHWS 


gwws 


h° 
H° 


cp a /2sw 


— Cf3 a /2cw 
S(3 a /2cw 


mwspa/sw 
m w cp a /sw 



Fig. 0.10: 



M = a g™*ggg*a B (gfA 1+9 ^A 2 )(B13) 



with the arguments for the integral functions as 
B(m 2 H± , nig., Trig.). The couplings are given in the fol- 
lowing table: 



Si Sj 




gzwSiSj 


9HSiSj 


h° H± 
H° H± 
h° G ± 
H° G ± 


Cf3a/2sw 
— Sp a /2sw 
Sf3 a /2sw 
C(3a/2SW 


— Cp a /2cw 
Sp a /2cw 

— Sf3 a /2cw 

— Cf3 a /2c W 


9H-h°H+ 

9h-h«h+ 

9H-h°G+ 

9h-h°g+ 



Fig. 0.11: 



ol gzHwsgzzs 



M = - ^ vv BoigfAt + g e v L A 2 ),(BU) 

s — mi, 



with the arguments for the integral functions as 
B(s,m s ,rn^). The couplings are given in the following 
table: 



s 


gzHws 


gzzs 


h° 


— Ci3a/2cw 

sp a /2cw 


m z s [Ja / s w c w 

m Z Cf3a/swCW 



Fig. 0.12: 

M = a 2 gwgwwsgHWS {[Co + {m s - u)D + (m^ - u)Di + (u - m 2 H± )D 3 ]A2 - 4:DiA 6 } , (B15) 

with the arguments for the integrals as C(m 2 , m 2 , s, m^, 0, m^), D(m 2 v ,ml,m 2 ,m 2 ^ ± ,u,s,m s ,rn 2 v ,0 7 rn 2 v ). The 
couplings are given in the following table: 



s 


gwws 


gnws 


h° 

H° 


mwspa/sw 
m w cp a /s w 


C{3 a /2SW 
— S[ia/2sw 



Fig. 



J w+ 



H+ (k 2 ) 



2vr 



gw [{mtVmb + m bymb) B i + ™ b y§ tb B ] , 



with the arguments for the integral functions as B(k , m?,m|), for k being the external momentum. 
Fig.|.2+l3: 



J w+ 



H+ (k 2 ) = 



[cpaSwgH- h°H+ (2-Bi (k 2 , m 2 h0 , m 2 H± ) + B (k 2 , to| , m 2 H± ) ) 



- spasw 9h- h° h+ {2Bi(k 2 , m 2 H a , m 2 H ±) + B (k 2 ,m 2 H a,m 2 H± )) 
+ s/ 3a s w g H - h o G +(2Bi(k 2 ,m 2 l0l m 2 v ) + B (k 2 , m%>, My/)) 



cp a Sw9H-H°G+ (2B 1 (k 2 ,m 2 H o,m 2 v ) + B (k? 



, ni H o . 



m w )) 



- cp a S(3 a m w (B Q (k 2 ,m\o,m%/) - Bi(k 2 , m 2 ha , m w )) 
+ Cf3 a Sf3 a m w (B (k 2 ,m 2 H0 ,mw) ~ B x {k 2 , m 2 H0 , m w ))] 



Fig. |.4: 



S G+if + ( fc2 ) = {{VmbVGtb + VHtbVGu)M™Z) + [mbmt(3/HfbVGib + VmbVctb) + <{VmbVltb + VmbVGtb)] A) 



N c a 



(B16) 



(B17) 



2tt 



+ ^(y^ttygtb + VmbVGtb)Bi} , 
with the arguments for the integral functions as B(k 2 ,m 2 ,m 2 ). 



(B18) 
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Fig. |.5+|6: 

S G+ff+( fc2 ) = £^ {9h- h°H+ 9h- h a G+ B o{k 2 \m 2 h0l m 2 H± ) + g H -H°H+gH-H°G+ B a{k 2 \m 2 H0 ,m 2 H± ) 
+ gH-h°G+9G-hOG+ B o{k 2 ,m^o,m^) + g H - h °g+9g- H°G+ B o(k 2 ,m 2 H0 ,m^) 



S -^^ [- (k 2 + m 2 )B (fc 2 , mlo , 4) + 2fc 2 5 1 (fc 2 , m 2 „ , m 2 w )] 



[(fc 2 + m 2 HO )B (k 2 ,m 2 H o,m 2 w ) ~ 1k 2 B x (k 2 , m 2 H0 , m 2 w )] j (B19) 



Fig. |B|.l: The real part of the W H mixing self-energy fixes the counterterms for the r yW + H and ZW + H 
vertices, following Eq. |lj: 

M = -^ T g v Re(J: w+H+ (m 2 H± ))(g e v R A 1 +g e v L A2), (B20) 
s — my 

where # 7 = -1, ,g z = s w /c w , and = ^w+h+ + = Rc ( 2 -w+ff+ ) + Hm(Y, w+H +). 

The W~ H~ counterterm cancels the real part of the W~H~ mixing self-energy, following the renormalization 
condition defined in Eq. [l^. Therefore, only the imaginary part of the W~H~ mixing self-energy contributes to 
Figs. [|.2 and [|.3. Similarly, only the imaginary part of the G~H~ mixing self-energy contributes to Fig. g.4; this is 
related to the imaginary part of the W~ H~ mixing self-energy by Eq. [l^. 

Fig. ||2: 

M = ^ a 9w n m (Z w+H+ (m 2 H ±))A 2 . (B21) 

m H± ~ m W 

Fig. |3: 



M = f^WW , _ m 2 )iIm{Ew+H+ {m 2 )){g eR Ai + geL Ai) (fi22) 



where Oyww = — 1 and gzww = —cw/sw- 
Fig. | .4: 



M = —. 4 ™ 9 T Gm h \ Hm(E w+H+ (m 2 H± ))(g^A 1 +g e v L A 2 ), (B23) 



where g lW G = ™>w and gzwc = -m z s w . 



APPENDIX C: MSSM CONTRIBUTIONS 

The matrix elements for e + e~ — > W + H~ in the full 
MSSM are given here for the first time. The sfermion tri- 
angle that enters Fig. 0.2 was computed in Ref. |40 , and 
agrees with our result. Fig. ||.l(b) with \i, X°j and xt m 
the loop is analogous to the top/bottom quark triangle 
diagram and can be checked by substituting top/bottom 
quark couplings for the gaugino couplings; the part in- 
volving the W coupling to left-handed gauginos is con- 
sistent with the top/bottom quark contribution given in 
the previous section. Formulae were given in Ref. [ ]38| for 
the W + H + and G + H + mixing diagrams, Fig. |[l-4; our 
results are consistent with theirs, although not enough 
detail was given in Ref. |3^] to check the overall signs of 
the diagrams. 

We define here our notation for the mixing matrices in 
the stop/sbottom and gaugino sectors and various cou- 



pling coefficients. 

The mixing in the stop sector is defined by: 



cosatti — sina t t2 
sbxottti + cos 




(CI) 



where th,R are the weak eigenstates and t\.2 are the mass 
eigenstates, and analogously for the sbottom sector. 
The chargino mass matrix is: 



X = 



M 2 \/2mw sin [5 
\/2rriw cos (3 /i 



(C2) 



which is diagonalizcd by the matrices U and V via 
VXW = M D . 

The neutralino mass matrix is: 
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Y = 



I Mi — mzsw cos P mzsw sin (3 ^ 

M 2 mzcw cos (3 —mzcw sin (3 

-mzsw cos (3 mzcw cos (3 — /i 

\ m^SvK sin (3 —mzcw sin/3 — /i / 



(C3) 



which is diagonalized by the matrix N via NY^N^ = No- fined as: 

The matrices that enter the W + XiXj couplings are 
defined as: 

0% = -l= Nii V^+N i2 V^ 0?j = -±=N* 3 Uj 2 +N* 2 U 3l . 



(C4) defined as: 



The matrices that enter the ZxJXj couplings are de- 



fined as: 



0. 



iL 



o. 



ill 



-VaV^--V i2 V* 2 + 6 ijS 2 w , 
-U* 1 U jl -\u* 2 U j2 + S ijS 2 w . (C5) 



O" L = -\N i sN* i + l -N ii N* A , Of = -Q'lf\ (C6) 



The matrices that enter the H XiXj couplings are 



1 



The matrices that enter the Zx\x^\ couplings are de- lowing table: 



Qij = A^/i + ^(N* 2 + N*i t<m6 w )V* 2 , 
Q'if = NisUj! - -j={N a + N a t a n6 w )U j2 . (C7) 

The if - //' coupling coefficients are given in the fol- 



J 



9He L u 
9Hd L u L 
9H-bJ L 
9H-b L i R 
9H-b a i L 



—raw sin 2(3j\/2sw 
— mw sin 

—mw sin2(3/^/2sw + (ml tan/3 + m 2 cot (3) j '^/2mw sw 
m t (fi + A t cot f3) j 'V2fflff sw 
mj(fi + Abtan f3) j '\J~2mw sw 
m t m b (ta.n (3 + cot (3) / \[2mwsw 



9H bl t 2 
dHb^h 
9Hb 2 t 2 



cosa t cosa b g H _i LiL + sinat cosa b g H _ bLiR + cosa t sma b g H _ bRiL + sma t sma b g H _ bRiR 
-sma t cosa b g H _- bLh , _ t ^^ Ui)H ^ bLtR 
-cosa t sma b g H -~ bLiL - sina t ^a b g H - bL i R 

sma t sma b g H - bLiL - cosat sma b g H _ bLiR - 



- cosa t cosa b g H - b j R - sina t sn\a b g H - bRt 



cosa t s\na b g H _ bRiR 
cos a t cos a b g H _ bRiL + sin a t cos a b g H _ bRiR 



■ sin a t cos a b g H - bRtL + cos a t cos a b g H - bR t R 



The G //' coupling coefficients are given in the following table: 



9Ge L C> 

9a-l L t L 
g G-b L i a 
9a-b R t L 



•9gM 2 

9Gb 2 t! 

9Gb ? M 



mw cos 2/3/^25^/ 
m^v cos 2/3 / V2sw 
mw cos 2/3 / 'v^siv + (m^ — m^/v^mwsw 
m t (— ii cot (3 + A t )/V2mwsw 
— m b (— [i tan /3 + / V^m^s^ 




cosa t cos a b g G -~ b j L +sina t cosa b g G _ bhlR + cosa t sma b g G _ bRlL 
-sina t cos a b g G _ b j L + cos a t cos a b g G - b j R - sinat sin a^.^ 
- cos a t sin a b .g G _^ - sina t sina b 5 G _^ fj? + cos a t cos a b g G _ bRlh 

sin a t sin a b g G - bi ij - cos a t sin ab9 G - br i R ~ sin a * cos a b9 G -b n t r , 



We now list our results for the MSSM diagrams. It is to be understood that diagrams involving charginos xt are 
summed over i = 1, 2 and diagrams involving neutralinos are summed over i = 1, . . . , 4. 
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Fig. |.l(a): 

M = ; 



{[(F-G + H)B + (m 2 Xk (F - G + H) - k x ■ k 2 H - m Xi m Xj m Xk J)C + (m 2 w F + ki ■ k 2 G - (h ■ k 2 - m^)H)Ci 
+ (-fci • k 2 F - m%± G-(fa-k 2 - m 2 H± )H)C 2 - 2(F + H)C 00 ] (gf'Ai + g e v L A 2 ) 
+ [HC + (-G + H)d + (F + 3H)C 2 + 2(F + H)(C 12 + C 22 )] [g e v R (A 3 + As) + 9v (A* + As)] 
+ [H vm C + (G im + + (F im + H im )C 2 ] (g e v R A 7 + g e v L A 8 )} (C8) 

where 

F = m Xi (g L H g^g w + g§9$gw), G = m Xj (g H g R g w + 9h9v9w), H = m Xk (g§g^g w + 9h9v9w), 
J = 9H9v9w+9k9v9w, 

F im = m Xi {9h9v9w ~ 9h9v9w)^ C = m X] {g H g R g w ~ 9h9v9w)> H im = m Xk {9h9v9w ~ 9h9v9w)- ( C9 ) 

The arguments for the integral functions are B(s, m Xj , m Xi ) , C (m w , s, m 2 H± , m Xk , m x . , m Xi ) . 
Fig. |^.l(b): Similar to Fig. |^.l(a), under the exchange of 



9h 9hj At, As terms flip sign. 
The couplings are given in the following table: 



(CIO) 





Xt X] Xk 




9* 


9 L z 


9 R z 


9 L h 


9 B h 


9w 


g W 


1.1(a) 
|l(b) 


Xi Xj Xk 

v9 

A.% Ay Aft 


-Sij 






Ofjswcw 
O'l^/swcw 


O'jf/swcw 
0'/ R /s w c w 


-Q'X* smP/s w 
-Q'i* sinP/sw 


-Q'kC cos P/s w 
-Q%* cos P/s w 


Of k /sw 


0*/s w 
Of k /s w 



Fig. 



M 



AN r a 2 



■9vs z s 3 9Hs i s h 9ws j s h {CooCffy^Ai + gy L A 2 ) 



(C 2 + C 12 + C 22 ) [g e v R (A 3 + As) + gf (A, + As)] } 



(Cll) 



with the arguments for the integral functions as C(m^, s, m 2 H± , m| , m|. , m|. ), N c = 1 for sleptons and N c = 3 for 
squarks. The couplings are given in the following table: 



Si 


s 3 


Sk 


gists. 


gzSiSj 


9HS.Sk 


gWSjS k 


V 


V 


e-L 


-gf 


-gf 


gHe L v 


-gw 


e-L 


e-L 


V 


gf 


gf 


gHe L i> 


gw 


ul 


UL 


d L 


-9? 


-gf 


9Hd L u L 


-gw 


d L 


d L 


UL 


gf 


gf 


9Hd L u L 


gw 


U 


tj 




-gfM^M^ - Q^MkMk 




9m k u 


~g w Ml 3 Ml k 


bi 


h 


tk 


gfM b Lt M b Ll +gfM b Rt M b Rl 


9z LM Li M Lj + 9z RM Ri M Rj 


gHbih 


g w Mt k M b Lj 



where M L1 = ML = cos at, — M L2 — M R1 = sina t , M b Ll — M R2 — cos a&, and — M b L2 — M R1 — sinaf,. 
Fig.|.3: 



M 



2 gw9 XieI g*t [(Co + C 2 )Tn x% g L H + C 2 m X] g R ] A 2 , 



(C12) 



with the arguments for the integral functions as C(0, it, rn 2 H± , m x . , m| , m Xj ) . The couplings are given in the following 
table: 



Xi Xj i 


g L H 


g§ 


g L f 


g R r 


X° x] e L 
xt Xj v 


~Q' R * sin ft a w 
-Q' :j R * sin P/sw 


-Q'i 1 -* cos p/sw 
-Q'ji* cos p/sw 


N^/V^sw + N^/Vlcw 
-Vil/aw 


-Uji/sw 
-N j2 /V2s w + N jl /y/2c w 
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Fig. 



M 



with the arguments for the integral functions as C(0, u, m H± , m~ L , m" , m~ 
table: 



U 2 gwg^ee L gxivC.9He L »C 2 A2, (CI 3) 

■ " J 1 ■ - - ■ ■ 1 -' 1 The couplings are given in the following 



Xi 


q L - 


q R - 


x° 


N? 2 /V2s w + Nfjy/2cw 


-N t2 /V2s w + Nn/VZcw 



Fig.|.5(a,c): 



M = a 2 g XzeI g XkeI {[m Xt g^g^(-C Q -ml ] D -m 2 w D 1 -u(D 2 + D 3 ) + 2D 

1 x,9w9h ("Co - ( u + ™ X] )D - (« + m^)Di - 2uL> 2 - (u + m 2 H± )D 3 + 2D 00 ) 



A 2 



+ m x 

+ m Xk9w9H (Co + m 2 Xj D Q + u{D x + D 2 ) + m 2 H± D 3 ^ + m Xt m Xj m Xk gwgh D o 

+ [™^9w9h (-2OD33 + ^13 + D 23 ) + m X] g w g% (-2D 3 - 2(D 33 + D 13 + D 23 ))] A 4 

+ [m Xt9w g L H (-2{D 2 + D 3 ) - 2{D 22 + D 33 + D 12 + D 13 + D 23 )) 

+ ™ X j9w9h (~2i?o - (2X>i + 4£ 2 + 4£ 3 ) - 2(£> 22 + £33 + D 12 + D 13 + 2D 23 )) 



with the arguments for the integral functions as C(0, 0, s, m 2 ^ , m~, m- Xi ), D(m 2 v , 0, 0, rn 2 H± , u, s, m^, , m Xk , wi? 
The couplings are given in the following table: 



(C14) 

™2 \ 



Xt Xj Xk i 



9x,el 



9 



Xkel 



9h 



9w 



9w 



0%/sw 

-Of k /aw 
-O^/sw 
-Oii/sw 



|.5(a) 
|5(b) 
&5(c) 
06(a) 
06(b) 



X° Xj + X° e L 

X? Xj + X° e fl 

xt X°j Xfe 5 

X° xt X°k e L 

Xi xf xl e R 



+ 



y2sw y2cw 

-V2Nn/c w 
-V? x /s w 



N k2 



+ 



v2sw y/2cw 

-V2Nn/c w 



y2sw y/2c\v 

-V2N* kl /c w 
—Vki/sw 

N k2 , N, 



+ 



~V2N* kl /c w 



-Q'if* smP/ Sw 
-Q\f* sinp/s w 
-Q'£* sin/3/sw 
-Q' k f sin (3/sw 
-Q'** sin/3/s w 



cos (3/s w 
-Q'^cos (3/s w 
-Q'fi*cos (3/s w 
-Q' k *cos (3/sw 
-Q' k L *cos (3/s w 



O^/sw 
O^/sw 
-Of k / Sw 

-Of} /aw 



-Op/aw 



Fig. ||.5(b): Similar to Fig. @.5(a,c), under the exchange of 

9h^9h> 9w^9w, A 2 <->Ai, A4 <-> A3, Ae <-> A 5 

Fig. 0.6(a): 



(C15) 



M 



°"\ J»* ,/{ 



m x,9w9H 



l Xi9v/9n 
.2 



(C 



D +m^ ±J D 1 +i( J D 2 + J D 3 ) 



^2 



(-Co - (m£. + t)D - (t + m 2 H± )D l - 2tD 2 - (t + m z w )D 3 + 2D 00 

+ m xk9w9n (-C - rn\.D - t{D x + D 2 ) - m^D?, + 2£> o) + m Xi m Xj m Xk g w 
+ [m Xi g w g L H 2D 3 + m Xj g w g^ {-2D - (4£>i + AD 2 + 2D 3 ) - 2{D n + D 22 + 2D 12 + D 13 + D 23 j) 
+ ™ Xk 9w9H (~2(X>i + D 2 ) - 2(D n + D 22 + 2D 12 + D 13 + D 23 ))] A 4 

+ [m X3 9w9H (-2£>i - 2(Dn + D 12 + D 13 )) + m Xk g w g%(-2)(D n + D 12 + D 13 )] A 6 } , (C16) 
with the arguments for the integral functions as C(0, 0, s, rn Xk , m? m 2 ,.), D(irig± , 0, 0, m^, t, s, m x . , m Xk , to? m Xi )- 



Xk'—l- 

Fig. |2|.6(b): Similar to Fig. g.6(a), under the exchange of 

9h^9h, 9w^9 W i A 2 <->Ai, A A ^A 3 , Ae <-> A 5 , 

Fig. 07: 

X = -2a 2 g^ eiL g R eSL g H s L 9gw[D 00 A 2 -(D 3 + D 33 + D 13 + D 23 )A4 
- (D 2 + D 3 + D 22 + D 33 + D12 + D 13 + 2D 23 )A 6 ] , 



(C17) 



(C18) 



with the arguments for the integral functions as D{m w , 0, 0, m 2 H ± , u, s, to~, Tn~ L , m x . , m 2 L ). The couplings are given 
in the following table: 
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Xi 


q L - 


g R - 






N t2 /V2s w +N ll /V2c w 



Fig. |.8: 



M = -2a 2 g^ e9 gK eCj g H e L i>gw[DooA2-(D 1 +D 2 + D 11 + D22 + 2D 12 + D 13 + D 23 )A4 
- (D 1 + D n + D a + D 13 )A 6 ] , 



(C19) 



where the arguments of the integral functions are D{m 2 ^± , 0, 0, mfy, t, s, mj^ , m|, m^,. , The couplings are given 
in the following table: 



Xr 






xi 


-V*Js w 


-Vn/s w 



Fig. 19: 



{Doo{m Xi g% + m x .g§)A 2 + [(D 3 + D 13 + D 23 )m x% g L H + (D 13 + D 23 )m Xj g%] A A 
+ [-(D 2 + D 12 + D 22 )m Xi g L H - (D 12 + D 22 )m Xjff g] yl 6 } , 



(C20) 



with the arguments for the integral functions as D{m 2 v , 0, m 2 H± , 0, u, t, m 2 , m-^ , m^, , m Xi )- The couplings are given 
in the following table: 



Xi Xj 


9 L h 


9§ 


q L - 






-Q'Jt sin p/sw 


-Q%* caa P/s w 


-Va/sw 


N ]2 /V2s w + Nji/V2cw 



Fig. |.10: 



M = a2 9x t eV9 Xj ee L 9He L 

{ m x t m Xj9w D o + 9w(~ c o- m Xj D o ~ m 2 w D x - uD 2 + 2D oa ) 
+ 2g*(D 13 + D 23 )A 4 - 2g*(D 2 + D 12 + D 22 )A} , 



A 2 



(C21) 



with the arguments for the integrals as (7(0, m^ ± ,t, m x .,m~, m\ L ), D(m 2 v ,0,m 2 ^ ± ,0,u,t,m Xj ,m 2 , i ,m 2 ,,rn 2 L ). The 
couplings are given in the following table: 



X* Xj 


5m/ 


9% 


„R . 

i>Xiev 


q L - 


■y~>~ -y® 
A 7 


-OfJs w 


-OfJsw 


-Vn/sw 


N* 2 /V2s w + N*jV2c w 



Fig. Ill: 



M 



N r a 2 



-gvws l s ] 9Hs,s ] B {g'y H Ai+ gyA 2 ) 



(C22) 



with the arguments for the integral functions as B(m 2 H± , m 2 s . , mg.), N c = 1 for sleptons and iV c = 3 for squarks. The 
couplings are given in the following table: 



Si Sj 


gHs l s 1 


V 




gHe L 9 


UL 


d L 


9Hd L u L 


h 


h 


gHbJx 


h 


b 2 


gm 2 t 1 


h 




guilt?. 




b 2 


9hI?m 



g^wSiSj 



(e e + e u )/V2sw 
+ e<i)/V2sw 
{e u + e d ) cos a t cos a b / V2sw 
- (e« + e d ) cos a t sin a b / V2s w 
~{e u + e d ) sina t cos a b /V2s w 
(e„ + e d ) sin a t sin a b j \/2s w 



gzws z s 1 



-(e e + e u )/V2c w 
-(e« + e d )/V2c w 
-(e u + e d ) cos a t cos ab/v^cw 
(e« + e d ) cos a t sin a b / \f2c w 
(e u + e d ) sina t cos a b / V2cw 
-(e„ + e rf ) sin a t sin ab/^cw 



Fig.|4|.l: 

s ^+H+( fc2 ) = ^{{^A9H9w+9§9w)+rn Xi (gjig^+g§g^)]B 1 +m Xi (g^+g§g^)B }, (C23) 

with the arguments for the integral functions as B(k 2 , m 2 ,., fn Xj ), for k being the external momentum. The couplings 
are given in the following table: 
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X* Xi 


9 L h 


9% 


9w 


9w 


Xi /Cj 


~Q'K* sin (3 /s w 


-Qij* cos (3 /sw 


0%/sw 


Ofi/sw 



Fig. |.2: 

Z b w+H Ak 2 ) = -^9HS i s j 9ws i s j ( 2B i + B o), (C24) 

with the arguments for the integral functions as B(k 2 , m|. , m|.), for fc being the external momentum. The couplings 
are given in the following table: 



Si Sj 


9HSiSj 


9ws l s 1 


V 




9He L u 


9w 


UL 


d L 


9Hd L u L 


gw 


tl 


h 




gw cos at cos a>b 


h 


b 2 




—gw cos at sin 


h 


h 


9Hbit 2 


— gw sina t cos ab 


*2 


b 2 


9Bb,t9. 


gw sin at sin a& 



Fig. |.3: 

S G+ff+( fc2 ) = {(9h9§ + 9H9h)M^ 2 Xj ) + h 2 x M9§ + 9Hgh)+^x,m Xj (g^gh + g^g§)] B 

+ k 2 (g L H g§ + gM)Bi}, (C25) 

with the arguments for the integral functions as B(k 2 ,m 2 , i ,m 2 , j ), for k being the external momentum. The couplings 
are given in the following table: 



Xi Xi 


9 L h 


9§ 


9 L a 


9g 


A,% A 7 


-Q'if* sin[3/s w 


-Qfi* cos /3/ s w 


-Q'tj sin P/sw 


Q[f cos /3/sw 



Fig. |.4: 

^ b G+H+ {k 2 ) = ^gHS&gas&Bo, (C26) 

with the arguments for the integral functions as B(k 2 , m|. , Wg.), for /c being the external momentum. The couplings 
are given in the following table: 



S{ Sj 


9HSiSj 


gaSiSj 


V 




9He L v 


gGe L i> 




d L 


9Hd L u L 


9cd L u L 


h 


h 


9m x t x 




h 




9Hb 2 t x 


9Gb 2 ti 


h 


bi 


9HbJ 2 


9GbJ 2 


t 2 


h 


9m 2 t 7 , 


9Gbii 2 
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